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ABSTRACT 

We present the first detailed photometric and spectroscopic study of the white dwarfs (WDs) in 
the field of the ~ 225 Myr old (logr c ; = 8.35) open cluster NGC 1039 (M34) as part of the ongoing 
Lick- Arizona White Dwarf Survey. Using wide-field UBV imaging, we photometrically select 44 WD 
candidates in this field. We spectroscopically identify 19 of these objects as WDs; 17 are hydrogen- 
atmosphere DA WDs, one is a helium-atmosphere DB WD, and one is a cool DC WD that exhibits 
no detectable absorption lines. We find an effective temperature (T e g) and surface gravity (log<?) 
for each DA WD by fitting Balmer-line profiles from model atmospheres to the observed spectra. 
WD evolutionary models are then invoked to derive masses and cooling times for each DA WD. 
Of the 17 DAs, five are at the approximate distance modulus of the cluster. Another WD with a 
distance modulus 0.45 mag brighter than that of the cluster could be a double-degenerate binary 
cluster member, but is more likely to be a field WD. We place the five single cluster member WDs 
in the empirical initial-final mass relation and find that three of them lie very close to the previously 
derived linear relation; two have WD masses significantly below the relation. These outliers may have 
experienced some sort of enhanced mass loss or binary evolution; however, it is quite possible that 
these WDs are simply interlopers from the field WD population. Eight of the 17 DA WDs show 
significant Ca II K absorption; comparison of the absorption strength with the WD distances suggests 
that the absorption is interstellar, though this cannot be confirmed with the current data. 
Subject headings: white dwarfs — open clusters and associations: individual (NGC 1039) 



1. INTRODUCTION 

This is the fourth installment in a series of papers in- 
vestigating the white dwarf (WD) content in young open 
clusters. The primary goal of these studies is to better 
constrain the integrated mass lost by low to intermedi- 
ate mass stars during their main sequence and advanced 
stages of evolution. This mass loss is most commonly 
represented as an initia l-final mass relation (IFMR, e.g., 
IWeidemannlll977l |2000D . Constraining the IFMR is cru- 
cial to understanding the chemical enrichment due to 
stell ar mass loss and sta r formation efficiency in galax- 
ies (jFerrario et al.ll2005l ). The shape of the IFMR has 
a direct effect on the shape of the WD luminosity func- 
tion (WDLF), which is commonly u sed to constrain the 
age of the Galacti c thin disk (e.g., IWinget et alJll987t 
lHarris et al.l |2006[). the thi ck disk, the Galactic halo 
' (e.g., 



_ von Hippel et al.ll2005bf). a nd globular clusters 
lHansen et aUl2002L I2004L 12007( 1 . 

Much work has been done to constrain the IFMR semi- 
empiricall y. The first purely photog raphic study was 
done by iRomanishin &; Angel ((19801 ): spectroscopy of 
WDs, which is necessary to determine their masses and 
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cooling times, wa s added to studies o f the IFMR in a 
series of papers bv lKoester fc Reimersl (e . g.,ll981lll996( ). 
Motivated by this progress, W eidemannl (|2000( ) updated 
his semi-e mpirical IFMR. M o re recent work on open 
clusters by[Claver et al. | 1(200 111. IKalirai et all (1200511 and 
IWilliams et "all (|2004af l led iFerrario et all ((2005( 1 to re- 
determine the empirical IFMR. They found that a lin- 
ear function adequately described the observational data. 
More recent purely empirical IFM R determinations (e.g., 
I Williamsll2007t IKalirai et al.ll2008f ) also find the empirical 
IFMR can be adequately parameterized as a linear func- 
tion; the slopes derived in these stud i es are consistent 
with the slope found by IFerrario et al . (2005). 

Another aim of these studies is to more accurately de- 
termine the upper limit on the main sequence mass of a 
star that evolves into a WD. This limit is an important 
parame ter in modeling dw arf galaxy evolution in par- 
ticular (|Dekel fc Silld Il986h . as it constrains the amount 
of energy and matter produced during the evolution of 
a given stellar populat ion. Recent theoretical studies 
(jPoelarends et al.l [2008) of the advanced stages of stel- 
lar evolution find that the lower limit on the mass of a 
solar metallicity star that will ignite carbon burning is 
~ 7Mq. Only stars with masses below this limit should 
be able to form CO-core WDs. Between ~ 7M ffi and 
~ 9M©, stars modeled in iPoelarends et all (j2008) ignite 
carbon but do not ignite neon and oxygen, so that some 
of them may form ONe-core WDs. These limits have yet 
to be tested observationally. WD studies to date have 
investig ated WDs with initia l masses between ~ 1 . 2 and 
6.5M (jFerrario et al.l 120051: IKalirai et al.ll2008f) . It is 
important to explore the WD populations of young open 
clusters with possible massive, still relatively bright WDs 
for this reason. 
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A third use of these studies of WDs in young 
open clusters is to compare cluster ages as deter- 
mined by main-sequence fitting with WD cooling times 
(|von Hippel et al. 2005a) . Consistency of these ages in- 
creases confidence in both stellar evolutionary models 
and WD cooling models. Consistency across a wide range 
of star cluster ages further increases confidence in studies 
seeki ng to derive star-for mation history from the WDLF 
(e.g jLiebert et al.ll2005h . 

Our general approach is to use photometric observa- 
tions of young clusters in order to select likely WD cluster 
member candidates based on their location in the color- 
magnitude and color-color (U-B versus B-V) diagrams 
for the fields. We then obtain high signal-to-noise spec- 
tra of these candidates. With spectroscopy of cluster WD 
candidates, it is possible to unambiguously determine if 
the objects are WDs. For the bona fide WDs, we deter- 
mine the effective temperature (T e g) and surface gravity 
(log g), and from those quantities and WD evolutionary 
models, we derive the cooling age (twd), mass, and lumi- 
nosity of each WD. For those WDs with cooling ages less 
than the cluster age and distance moduli consistent with 
cluster membership, subtraction of the WD cooling age 
from the age of the open cluster results in the lifetime of 
the progenitor star, and stellar evolutionary models can 
then be used to determine the progenitor star mass. 

As part of our ongoing program to identify and spec- 
troscopically analyze the WD populations of open clus- 
ters, the Lick- Arizon a White Dwarf Survey (LAWDS, 
iWilliams et"aT][2004ah . we have obtained imaging data 
to identify candidate WDs in NGC 1039 and high signal- 
to-noise spectra of WDs in the field of NGC 1039. 

In this paper, we ass ume the solar metallicity value 
of lAsplund et al] (I2004L Z = 0.0 13) and the extinction 
curve of iRieke fc Lebofskvl (|1985f) with R v = 3.1. 

1.1. Previous Studies of NGC 1039 

The fi rst study of WD s in N GC 1039 (M34) was 
done bv lAnthonv-Twarod (|1982D . In this study, UBV 
photographic plates were used to obtain photometry 
for stars in the cluster field and to select objects with 
significant blue color excess. The color-magnitud e se- 
quences for WDs de rived by iSion &: Liebert] (|1977l ) and 
iKoester et al.l ()1979l ) were compared with these observa- 
tions. Four objects were found to be close to the WD 
cooling curves used in color-magnitude and color-color 
space, and so could be interpreted as being cluster mem- 
ber WDs. The number of field WDs of different apparent 
V magnitudes expected to appear in the ph otometry wa s 
calculated using the then-current WDLF (|Greenl ll980D . 
This number of expected field WDs (~ 23 with V < 21) 
was, in fact, larger than the number of observed blue 
objects (11 with B — V < 0.4), leading to the conclu- 
sion that NGC 1039 did not have a significant number of 
WDs. 

A number of recent studies h ave been done on NGC 
1039. Ilanna &T Schlcmmcr (1993) present proper motions 
to determine which objects are cluster members. They 
determine (m — M)o = 8. 28 and an age of log T r .i = 8.40 
using E(B - V) = 0.07. iMevnet et all (fl993D compare 
new stellar models to previously obtained data. They use 
E(B - V) = 0.10 an d calculate (m - M) n = 8.34 and an 
age of log t c i — 8.25. Uones fc Prosserl ( |1996l ). again using 
a proper motion selected sample, use E(B — V) — 0.07 



and obtain (m - M) = 8.38 and logr c/ = 8.30 - 8.40. 

The most re c ent s tudy of NGC 1039 is that of 
ISaraiedini et al] (|2004D . These authors combine UBVRI 
photometry from the WIYN open cluster survey with 
JHK photometry from 2MASS to obtain accurate dis- 
tance estimates. They assume E(B — V) = 0.10 and 
[Fe/H] = +0.07 and derive (m - M) = 8.67 ± 0.07. 
The spectroscopic metal abundance of N GC 1039 has 
been d etermined to be slightly super-solar. ISchuler et al.l 
(2003) find [Fe/H] = +0.07 ± 0.04, and new work from 
the WIYN Open Cluster Study finds [Fe/H] « +0.08 (A. 
Steinhauer, private communication, 2008). In this paper, 
we assume solar metallicity for the cluster. 

2. PHOTOMETRIC OBSERVATIONS & ANALYSIS 
2.1. Photometric Measurements and Calibration 

UBV imaging of a field centered on NGC 1039, as well 
as of two fields 25' to the north and south of the cluster 
center, was obtained in 2004 September and 2005 Decem- 
ber with the KPNO 4m MOSAIC camera. Exposures of 
5, 30, and 3 x 120 s were obtained at the northern field 
in both B and V and at the central field in V; exposures 
of 5, 30, and 3 x 150 s were taken at the central field 
in B. We obtained exposures of 5 and 120 s in B and 
V at the southern field. Exposures of 5, 30, 120, and 
3 x 900 s were obtained in U at the central field; all of 
these but the 30 s exposure were obtained in the northern 
field. The southern field was observed with exposures of 
15, 120, and 3 x 600 s in U. The northern and south- 
ern fields are positioned at very nearly the same right 
ascension (R.A.) as the central field. Each flanking field 
overlaps the central field by ~ 11' on either its north- 
ern or southern end. Seeing ranged between 078 and l'.'O 
in the F-band in 2004 September and was ~ 1'.'7 in the 
B-band in 2005 December. 

We reduced the data using the IRAF mscred mo- 
saic reduction package. The data were bias subtracted, 
trimmed, flat-fielded and projected onto the tangent 
plane. Object detection and aperture photome try were 
perfo rmed using the DAOPHOT II program (jStetsonl 
1987); aperture corrections were determined an d applied 
using the program D AO GROW (|StetsoiJll990h . 

Star-galaxy separation was performed by comparing 
the y-band flux in two successively larger apertures for 
each object: an inner aperture with a radius equal to the 
full-width half-maximum (FWHM) of the t/-band point 
spread function (w 0'.' 8) and an outer aperture with a 
radius of twice the FWHM. Stars were observed to form 
a tight sequence with a constant difference in aperture 
magnitudes of rs 0.3 mag for V < 23; we therefore iden- 
tified objects with aperture magnitude differences < 0.5 
and V < 23 as point sources. We consider only these 
point sources in the subsequent analysis. 

As the MOSAIC images of NGC 1039 were taken under 
non-photometric conditions, we calibrated these images 
using previously obtained UBV imaging of the central 
portion of NGC 1039 taken with the Nickel 1-m tele- 
scope at Lick Observatory on 2002 July 6. These images 
were reduced and aperture photometry was determined 
for each star. The Nickel imaging was calibrated using 
standard stars from iLandoltl (|1992h . The stars with cal- 
ibrated Nickel photometry were then used as secondary 
standard stars for calibrating the MOSAIC imaging. 
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Calibrations were obtained by solving the following 
transformation equations for both the MOSAIC and 
Nickel imaging: 

u = U+2.51o g t cxp +A a +A 1 (X-1.25)+A 2 (U-B)+A 3 (U-B) 2 (1) 
6 = B+2.5logt^ p +B +B 1 (X-1.25) + B 2 (B-V) (2) 
v = V+2.5logt eKp +Co+Ci(X-1.25)+Ca(B-V), (3) 

where u, b, v are total instrumental magnitudes; U, B, 
V are standard magnitudes; t exp is the exposure time; 
and X is the airmass. 

Zero points and color terms were calculated for the 
MOSAIC data; as local secondary standards were used, 
the airmass terms were not necessary. Values for the 
transformation coefficients are given in Table [1] The B 
and V color terms are consistent with average values de- 
termined bv lMassev et all (|2006l ) for the MOSAIC cam- 
era; the [/-band color terms are not directly comparable 
because we have chosen a quadratic formula. 

We compared our c alibrated MOSAIC phot ometry of 
46 objects with that of iJones fc Prosserl (|1996h . There is 
no significant trend in the magnitude of B — V color offset 
with magnitude or color between the two systems. The 
mean offset in V magnitude is 0.021, with a dispersion 
of 0.048. The mean offset in B~V color is 0.008, with a 
dispersion of 0.050. 

The left most panels of Figures [1] and [2] show our fi- 
nal calibrated photometry for the cluster in B — V, V 
and U — V, V color-magnitude diagrams. Note that 
sources with V < 13 are saturated in our photometry. 
IJones fc Prosserl ((1996) photometry is overplotted in the 
second left panel in Figure [T] for stars with a high prob- 
ability of cluster membership (> 67%, where probability 
is determined from their proper motion study). 

2.2. WD Candidate Selection 

We select WD candidates using photometric cri- 
teria based on synthetic evolutionary sequences of 
WDs kindly provided by P. Bergeron 5 ; these se- 
quences repres e nt an extension of those published in 
iBergeron et al.l (119951). including improved color calcu- 
lations (jHolberg fc Bergeron! 120061) and the more rece nt 
WD evolutionary calculations of Fontaine et al. (2001). 

The second right panel of Figure [1] and the mid- 
dle panel of Figure [2] show isochrones and WD cool- 
ing curves adjusted to the cluster's (m — M)y and 
E(B — V) using cluster parameters from IJones fc Prosserl 
(1996). The right-hand panels of both figures show 
the same isochrones and WD cooling curves, adjuste d 
using cluster parameters from ISarai edini et all (2004 ). 
The solid isochrones are Padova models for logr c / = 
8.35 ([Girardi et al.ll2002f) . interpolated for Z = 0.013. 
The dashed isochrones are Yonsei-Yale logr c ; = 8.35, 
Z = 0.013 models interpolated from solar-scaled abun- 
dance models using the code p rovided by the Yonsei-Yale 
Group (jDemarque et al.ll2004h . The dotted isochrone in 
Figure [T] is the fiducial main sequenc e of the sub-solar 
metallicity cluster NGC 2168 given in [von Hippe fet al.l 
([20021) . corrected to the metallicity of NGC 1039 using 
the method described in lPinsonneault et aTl (fl998). The 
cooling curves for DA and DB WDs of different masses 
are also shown. 

5 http:/ /www. astro. umontreal.ca/~bergeron/CoolingModels/ 



It is evident from Figure [T] that our photometry 
matches the photometry in the literature acceptably well, 
and that modeled isochrones follow the B — V versus V 
main sequence in our data. The best fit of the ob- 
served main-sequence to model isochrones is obtained 
when we adopt the E(B — V) and (m — M)y values from 
ISaraiedini et al.l (|2004l ); model isochro nes shifted to the 
distan ce modulus and reddening from IJones fc Prosserl 
(1996) fail to match the observed main sequences. We 
therefore simply ad opt the values o f E(B — V) and 
(m — M)v given in ISaraiedini et al.1 (|2004fl for the re- 
mainder of our analysis. 

As seen in Figure [21 not one of the plotted isochrones 
appears to match our main sequence photometry in the 
fainter parts of the U — V versus V color-magnitude 
diagram well; the observed main sequence is signifi- 
cantly redder than any theoretical isochrone for V > 14 
(My > 5). The cause of this discrepancy is unclear. 
Given the agreement between models and observations in 
the B — V,V plane, the problem is likely either shortcom- 
ings in [/-band model spectra for cool stars or calibration 
issues in our [/-band photometry. Given the agreement 
between our WD photometry and photometric models 
(see N4.2p , the latter seems unlikely. Figure [3] shows a 
color-color diagram with our photometry in the field of 
the cluster. 

Figures H] and [5] show the sections of the color- 
magnitude diagrams and color-color plot surrounding 
the WD cooling sequences with our photometry of the 
field of the cluster. Cooling curves are shown as in Fig- 
ures [TJ [21 and [3] The boundaries of our selection re- 
gions for cluster WD candidates are shown with dashed 
lines. The ^-magnitude limits for WD candidate selec- 
tion are 13.69 < V < 21.38. We also required that all 
WD candidates have B-V < 0.30 and U - B < -0.41. 
These boundaries include all photometric models for DA 
WDs in the mass range 0.4M Q < M < 1.2M with 
twd ^ 8.48, adjusted for the adopted distance and red- 
dening of the cluster. Objects which meet all of these cri- 
teria were selected as possible WD candidates for spec- 
troscopic followup. They are marked in Figures 2] and 
[5] with open circles. Table [2] lists UBV photometry 
of all 44 selected WD c andidates. Seven objects from 
lAnthony-Twarogl (|1982f ) are recovered. See our §5.1.11 
for a more detailed description of these objects. 

3. SPECTROSCOPIC OBSERVATIONS & ANALYSIS 
3.1. Spectroscopy of Candidate WDs 

Spectroscopic observations of selected WD candidates 
were obtained over several observing runs between UT 
2001 August 22 and 2007 Jan uary 20 with the bl ue cam- 
era of the LRIS spectro graph ([Cohen et al.lll994h on the 
Keck 1 10-m telescope (|Oke et al.lll995l) . The 2001 Au- 
gust observations used the initial LRIS-B engineering- 
grade 2048 x 2048 SITe CCD; all other observations used 
the New Blue Camera, consisting of two 2k x 4k Mar- 
coni CCDs. The 2001 observations therefore have lower 
sensitivity in the blue. 

A l"-wide longslit at parallactic angle was used with 
the 400 1 mm" "\ 3400 A blaze grism for a resulting spec- 
troscopic resolution of ~ 6 A. The spectra were reduced 
using the onedspec package in IRAF. Overscan regions 
were used to subtract the amplifier bias, and cosmic rays 
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were removed from the two-dimensional spectra using 
the L.A. Cosmic La placian cosmic ray rejection routine 
(|van Dokkuml l2001). We then co-added any multiple ex- 
posures of an individual object and extracted the one- 
dimensional spectrum. We applied a wavelength solution 
derived from HgCdZn arc lamp spectra. We applied a 
relative flux calibration determined from longslit spectra 
of spcctrophotomctric standard stars. 

Spectra of WDs are shown in Figure [5] Spectroscopic 
identification of each WD candidate is given in Table 
The major non-WD contaminants in the sample are ac- 
tive galactic nuclei (AGNs). To determine redshifts for 
the latter objects, we cross-correlated the spectra with 
the Sloan Digital Sky Survey (SPS S) composite quasi- 
stellar object (QSO) spectrum of IVanden Berk et al.l 
(|2001f) as described in I Williams et all (|2004bl ). 

3.2. WD Parameter Determination 

T e g and logg were determine d for each WD usin g 
simultaneous Balmer-lin e fitting dBergeron et alJll992f ). 
applied as described in IWilliams fc Boltel (|2007f ) . The 
model spectra used for the fits are from Koester's model 
grids; d etails of the input p hysics and methods can be 
found in lKoester et al.l (|200lD and references therein. DA 
evolutionary models were pro yided freely by P. Berg- 
eron, and con sist of models o f IWoodl (|1995h (for T c g > 
30,000K) and iFontaine et all (|2001l ) ffor T eff < 30,000 
K), with "thick" hydrogen layers (qn = 10~ 4 ) and mixed 
C/O composition. These models were used to calculate 
the mass (Mwd) and cooling age (twd) of each WD. 

We consider three primary sources of error in these 
parameters: internal fitting errors (what is the expected 
distribution in the fit values if we fit multiple observa- 
tions of the same star with the same signal to noise?), 
external fitting errors (what is the expected distribution 
in the fit values if others were to observe the same star 
and derive independent fits?), and systematic error due 
to the cluster age, which affects all WDs in the same 
sense, and so should not be considered as a random er- 
ror to be added in quadrature. We discuss the two fitting 
errors here, and the systematic errors in §4.21 and §5.41 

Internal fitting errors were determined empirically. 
The noise measured for each spectrum was added to the 
best-fitting model spectrum convolved with the instru- 
mental response. These simulated spectra were then fit- 
ted by the same method; nine iterations were used to 
calculate the scatter in T g, log g, Mwd, and twd- The 
errors thus determined (<7T off ,int and cri og g,int) represent 
only random errors due to observations and spectral fit- 
ting. 

External fitting errors can be estimated for our data 
fro m the comparison of field WD parameters performed 
in IWilliams fe Boltel (|2007l ). This analysis found that 
the scatter in WD parameters due to the use of different 
instrumentation, spectral fitting routines, and model at- 
mospheres is er(T c ff) re HOOK and cr(logg) re 0.12. These 
are added in quadrature with the internal fitting errors 
to produce the adopted values for <7T cff , tot and criog g ,tot 
for each WD. 

The atmospheric fits and derived WD masses and ages 
are given in Table [H The Balmer-line fits are shown in 
Figure [7] A distance modulus for each WD was derived 
by comparing the observed V magnitude to the absolute 
magnitude My calculated from the best-fitting model 



atmosphere and the appropriate WD cooling model. 

4. WD CLUSTER MEMBER IDENTIFICATION 

4.1. WD Cluster Membership Determination 

We consider two criteria when determining whether 
a WD is a potential member of NGC 1039. (1) Its 
(m — M)v must be between the tw o determinations 
of the clus ter (to — M)y made by iJones fc Prosserl 
(11996[ ) and ISarajedini et al.l (|2004h and (2) its cooling 
age must be less than the cluster age. Figure [8] shows 
(to — M)y for each observ ed D A WD in a hist o gram . 
The IJones fc Prosserl (|1996Q and lSaraiedini et alj (|2004f ) 
distance moduli determinations are marked with dashed 
lines. We select five DA WDs to be potential members 
of NGC 1039 based on this distance modulus and age 
selection. Those WDs with dist ance moduli up to 0.75 
V magnitudes brighter than the IJones fc Prosserl (|1996l ) 
cluster distance modulus are in the appropriate bright- 
ness range to be in double-degenerate binary systems in 
the cluster. There is one such candidate among our DA 
WDs which is younger than the cluster. Table Q] lists all 
selected cluster member candidates. 

4.2. Progenitor Mass Calculation 

The progenitor mass (initial mass) for each candidate 
cluster member WD was calculated by first subtracting 
the WD cooling age (twd) from the cluster age. The 
resulting age difference is the total lifetime of the pro- 
genitor star from the zero-age main sequence (ZAMS) 
to the thermally pulsing asymptotic giant branch (AGB) 
phase; the amount of time between the AGB phase and 
the entrance of the star into the WD cooling track is as- 
sumed to be negligible. We then calculated the ZAMS 
mass of a star with that total lifetime fr om the solar- 
met al licity stellar evoluti o nary models of iGirardi et alJ 
(20001 and lBertelli et al.1 (fl994). The progenitor mass 
for each candidate cluster member DA WD is given in 
Table [4] for three ages of NGC 1039 spanning the range 
of recently published ages: logr c ; = 8.30 (re 200 Myr), 
logr ; = 8.35 (re 225 MyjQ , and log r c/ = 8.40 (re 250 
Myr; IJones fc Prosser| [l996) . Errors in the initial masses 
are la variations resulting from the combined internal 
and external fitting errors in log g and T c g. 

Figures and [TO] show our WD candidate selection 
region of color-magnitude and color-color space with ob- 
jects spectroscopically identified. Of particular note are 
the filled and open circles. Large filled circles repre- 
sent high mass (> O.8M ) DA WD cluster members 
(LAWDS 15, 17, and S2; see the next section). These lie 
quite close to the plotted O.8M DA WD cooling curves 
in both color-magnitude diagrams in Figure [5] Two of 
these lie close to the plotted cooling curves in the color- 
color diagram (Figure IT0|): one of them (LAWDS 17) lies 
to the right of the cooling curves. The intermediate- 
sized open circle represents the object selected as a can- 
didate double-degenerate cluster member (LAWDS SI). 
Intermediate-sized filled circles represent lower mass DA 
WDs selected to be cluster members (LAWDS 9 and 20). 

4.3. Field WD Contamination 

We select WDs for cluster membership based only on 
their distance and cooling time, yet it is quite possible 
for field WDs to satisfy these criteria as well. We esti- 
mate the probability of finding field WDs meeting our 
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selection criteria in the WD selection volume using the 
lHarris et~a l. (2006) luminosity function. In the 36' x 86' 
MOSAIC field, ~ 2.4 field WDs with M v < 11.8 (the 
faintest My of the confirmed WDs) are expected. 

Three of our WDs selected as cluster members are 
~ 0.8M©, roughly the value that would be expected for 
NGC 1039 WDs based on the cluster age and previously 
published empirical IFMRs. (Another WD of similar 
mass, LAWDS S3, is too close [(to - M) v = 7.72] to 
be associated with the cluster.) As ~ 22% of field W Ds 
are more massive than O.8M ()Liebert et al.l [20051 ). we 
would only expect to find an average of ~ 0.5 massive 
WDs in a random field, as opposed to the three we have 
identified. The probability of finding three high mass 
field WDs in our field is small (P ~ 0.015, assuming 
Poisson statistics); it therefore seems likely that these 
WDs are bona fide cluster members. 

Two of the five possible cluster members (LAWDS 
9 and LA WDS 20) have masses typical of field WDs 
(~ O.56M0 iBergeron et al.l [1992). indicating that these 
may be interlopers from the field; two interlopers would 
be completely consistent with the expected number of 
field WDs calculated above. It therefore seems plausi- 
ble that both lower mass candidate cluster WDs are not 
true cluster members, but we cannot be certain of this 
with our current data. Proper motion measurements are 
needed to clarify the membership status of these WDs. 

However, as we are attempting to construct a purely 
empirical IFMR, it is dangerous to throw out LAWDS 
9 and LAWDS 20 simply because they do not meet our 
preconceived notions of what the IFMR ought to be. We 
now examine some possible explanations for how these 
two WDs can be cluster members and yet have relatively 
low masses. 

1. The WDs may be members of binary systems and 
may have a lowered final mass due to mass transfer 
onto their binary partner. This scenario has been 
invoked to explain the lowest mass WD s in the field 
WD population (e.g.. iKilic et al.ll2007h . Although 
these potential NGC 1039 WDs are significantly 
higher mass, a similar mechanism could still oper- 
ate. However, as the colors of these two WDs are 
not abnormal for single WDs, the companion re- 
sponsible for the binary evolution either must be 
a very faint compact object, a very low mass star 
or brown dwarf, or the companion must have been 
lost. 

2. iKalirai et alj (|2007f ) have shown via spectroscopic 
measurements that many WDs in the old, super 
metal rich open cluster NGC 6791 are low-mass 
helium-core WDs, and have posited that the high 
metallicity of the star cluster may have led to en- 
hanced mass loss during the post-main-sequence 
evolutionary phases of these stars. A similar mech- 
anism could be responsible for the apparently low 
masses of LAWDS 9 and LAWDS 20. 

However, the metallicity of NGC 1039 is consistent 
with solar (or slightly super-solar) metallicity. Fur- 
thermore, its metallicity is significantly lower than 
that of the Hyades and Praesepe clusters, and no 
abnormally low mass WDs are seen in Hyades or 
Praesepe. Enhanced mass loss therefore does not 



seem a likely explanation for the low mass of these 
two WDs. 

Given these points, we consider the most plausible ex- 
planation for the comparatively low masses of LAWDS 
9 and LAWDS 20 to be that these are interlopers from 
the field WD population. However, without proper mo- 
tion membership determinations, we cannot rule out that 
these stars are indeed cluster members. 

5. DISCUSSION 
5.1. Notes on Individual Objects 

Out of 44 photometrically selected WD candidates, we 
spectroscopically identified 32 objects. Twelve objects 
are AGNs, one is an A- type star, and 19 are WDs. Sev- 
enteen of these WDs are hydrogen-atmosphere DAs; of 
these, five are selected to be cluster member DA WDs, 
and an additional DA WD may possibly be a double- 
degenerate cluster member. 

LAWDS 26: This WD exhibits He I absorption lines in 
its spectrum and no hydrogen lines. We therefore classify 
it as a DB WD. The equivalent width of the He I 4471 
A absorption l ine is 13.4 A, suggesting T e s ~ 15, 000K 
for log g = 8.0 (Kocstcr 1980). A detailed analysis of this 
object and its probability of cluster membership will be 
presented in a future paper. 

LAWDS 7: This object shows no sign of any atmo- 
spheric features in the observed spectral range, so we 
classify it as a DC WD. As such, it is certainly too cool 
to have been born during the lifetime of NGC 1039, so 
we do not consider it a cluster member. 

LAWDS SI: Double degenerates, being a potential 
source of Type la supernovae, are of particular interest. 
At low spectral resolutions, double DA WDs appear to 
be single objects with T c g a nd log q intermediate t o those 
of the two component WDs (IBergeron et al.|[l990f) . Clus- 
ter double degenerates would also appear overluminous 
(to a maximum of w 0.75 mag), and would be identified 
as foreground objects according to our membership cri- 
teria. It is possible that the DA WD LAWDS SI may 
be a double degenerate. It is ~ 0.1 to ~ 0.5 mag over- 
luminous if at the cluster distance. A lower WD mass 
could result from close binary evolution; its initial mass, 
assuming that it is a cluster member, is given in Table 
[4j As the mass of this object (0.58M©) is near the peak 
of the field WD mass distribution, it seems more likely 
that this WD is just a solitary field WD. 

LAWDS S2: We observe a very weak (EWw 0.5 A) 
He I absorption line in the spectrum of the DA LAWDS 
S2 at 4471 A. No other He I or He II lines are observed. 
The temperature of this WD is near the red edge of 
the so-called "DB gap," the temperature region between 
28,000 K and 45,000 K where very few helium-rich atmo- 
sphere WDs are known to exist. One explanation for the 
DB gap is that, in this temperature range, small resid- 
ual amounts of hydrogen in a helium-atmosphere WD 
form a thin veneer obscuring the helium atmosphere; hy- 
drogen layer masses as small as 10~ 14 M© are sufficient. 
At 28,000 K, convection reaches the surface layers, mix- 
ing the hydrogen back into the dominant helium layer 
(jFontaine fc Wesemaell fl987h . LAWDS S2 may, as in- 
dicated by its temperature being near the cool end of 
the DB gap, be starting its transition from a DA back 
into a DB spectral-type object. Alternatively, LAWDS 
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S2 m ay be a binary WD, w ith DA and DB components 
(e.g.. IWesemael et al.lll994l ). Additional spectroscopy is 
needed to permit further analysis of this potentially in- 
teresting object. 

5.1.1. Objects Recovered from [Anthonv-Twarod (198li) 

In the early 1980s, studies of WD s in open clusters were 
pione e red by several autho r s (e.g ., iRomanish in fe Angel 
Il980t iKoester fc Reimers! 119811) . lAnthonv-Twarog 
(fl98l included NGC 1039 in her study. Of the four 
ob jects identified as p ossible WD members of NGC 1039 
bv lAnthonv-Twarod (|1982ft . only one, A35149 (LAWDS 
30 in this work), meets our WD candidate selection cri- 
teria. It has been spectroscopically identified as an AGN 
with a redshift z = 0.711 ± 0.003 assigned from a sin- 
gle emission line (Mg II) identification . The other three 
candid ate WD cluster members from I Anthony- Twaro"g| 
(1982) are recovered in our photometry but are too red 
to meet our WD candidate selection criteria. We do se- 
lec t several objects as candi date WDs that were identified 
bv l Anthony- Twarogj (|1982h but not selected in that work 
as candidate WDs. These objects are A43036 (LAWDS 
33, not observed spectroscopically), A42085 (LAWDS 41, 
an A-type star), A16073 (LAWDS 26, a DB), A22189 
(LAWDS 22, a DA), A24144 (LAWDS 20, also LB 3570, 
a DA), and A15118 (LAWDS 9, also LB 3567, a DA). 

5.1.2. WDs with Ca II K Absorption 

Eight of the 17 DA WDs (47%) show significant Ca II 
K absorption, indicating that these are potential DAZ 
WDs. Measured equivalent widths (EWs) of the Ca II 
K absorption lin e s for all WDs are given in Table \5[ 
IZuckerman et al.1 (|2003h find that ~ 25% of DA WDs 
show intrinsic Ca II absorption; our rate of incidence is 
significantly higher than this. In addition, our EWs are 
sig nificantly larger th a n thos e observed in the majority of 
the IZuckerman et al] (|2003[ ) sample; the higher average 
temperature of our WD sample would imply enormous 
abundances for a large fraction of WDs, which seems im- 
plausible. 

There is a distinct possibility that the detected Ca II 
is interstellar in origin; our spectral resolution is too low 
to resolve individual clouds in the interstellar medium 
(ISM) or to measure any velocity differences between the 
Balmer lines and the Ca II K line. While there appears 
to be no significant correlation between the Ca II K EW 
and the derived WD distance modulus, the measured 
EWs fall within the scatter of the inters tellar absorp- 
tion m easured in the spectra of hot stars bv lMegier et al.l 
(2005, see Figure [TTj) . However, the candidate DAZs do 
not seem to show any projected spatial correlations, and 
there is no correlation between the position of these WDs 
and WDs that do not exhibit Ca II K absorption in the 
field wit h areas of enhance d diffuse Ha emission from the 
maps of lFinkbeinerl (|2003h . 

In short, we cannot determine conclusively whether the 
observed Ca II K absorption in the eight potential DAZ 
WD spectra is intrinsic to the WDs, or due perhaps to 
patchy interstellar clouds, though the interstellar origin 
seems most likely. Higher-resolution spectroscopy is nec- 
essary to further investigate the potential metal content 
of these WDs. 

5.2. Expected Number of Cluster WDs 



The number of expected cluster WDs may be calcu- 
lated by assuming an initial mass function (IMF) for the 
cluster and an upper limit on the mass of st ars which 
form WDs; such calculations are detailed in IWilliamsl 
(2004). We calculate the expe cted WD popula t ion o f 
NGC 1039 using this method. Uones fc Prosseri (|1996l ) 
combine the results from their proper motion study of 
NGC 1039 with those of llanna fc Schlemmer! (fl993T) to 
construct a luminosity function for the central ~ 4 par- 
sees of the cluster. They find 71 cluster members with 
My < 4.0, the completeness limit of their study. We 
assume a Salpeter IMF and an upper mass limit on WD 
progenitors of 8Af©, and find that we expect the cluster 
to contain 11.6 WDs in this central region. Assuming a 
binary fraction of 0.5 and a random distribution of bi- 
nary mass ratios, 5.5 WDs should be detectable i n the 
~ 4 parsec field imaged by Uones fc Prosseri (1996); the 
r est will be hidden i n bin ary stars. 

Uones fc Prosseri ()1996l ) assume a projected surface 
density profile for the cluster of the form p(r) = poe^ r ^ r °\ 
where p(r) is in units of stars per square arcminute. They 
calculate r , the scale length for the projected surface 
density of cluster members, for several different V mag- 
nitude ranges. We normalize this projected surface den- 
sity profile usi ng the number of WD s expected in the 
area imag ed by J ones fc Pro sscr (1996, calculated above) 
and integrate it over the area of our MOSAIC imag- 
ing. We choose the largest value of ro = 9.'1 in order 
to avoid underestimating the expected number of WDs 
in the cluster. The result of this integration, 7.7 WDs, 
is the expected number of cluster member WDs in the 
field imaged in this study. If all six WDs observed at the 
cluster (m — M)v are cluster members, this result is con- 
sistent with our findings. If only the three massive WDs 
are cluster members, then NGC 1039 may be deficient 
in WDs; this deficie ncy has been ob served in other poor 
open clusters fe.g.. IWilliamsl 120041 ) . We also integrate 
the projected surface density profile to a radius of infin- 
ity and find that the absolute maximum number of WDs 
expected is 9.9. For smaller values of ro, the expected 
number of WDs decreases. 

5.3. Cluster Distance and Reddening Determination 
from WDs 

As noted previously, there is considerable variation in 
values of the distance modulus of NGC 1039 reported in 
the literature. An alternative method of estimating clus- 
ter distances involves measuring the distances of WDs as- 
sumed to be cluster members. The results of this method 
are independent of any assumptions we have made about 
the distance of the cluster, other than the membership 
distance criterion (which involved a broad selection). 

We first limit the analysis to the three massive 
(Mwd > 0.8M Q ) cluster members; the weighted av- 
erage of (m — M)v of these objects is 8.916 ± 0.133. 
This value matches the Saraiedi ni et al.l (|2004l ) determi- 
nation of (m — M)v = 8.98 ± 0.06 within the error bars. 
If we include in the calculation the lower mass poten- 
tial cluster members, we obtain a weighted average of 
(m — M)y = 8. 8 66 ± 0.100, also in agreement with the 
Saraiedi ni et al.l (|2004l) v alue. 



Saraiedini et al.l (|2004l ) use near-IR photometry from 



2MASS in addition to optical photometry in or- 
der to determine a distance modulus, as opposed to 
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Pones k, Prosserl (|1996|), whose study i s confin ed to BVI 
optical photometry. ISaraiedini et al.1 (|2004f ) note that 
when the optical phot ometry alone is fitte d by mod- 
ern isochrones (as in iRaffauf et al.l (120011) . who use 
iGirardi et all (|2000l ) and|YLe±al| (j2001|) isochrones), a 



distance m odulus consisten t with their measurement is 
obtained. ISaraiedini et al.l (|2004f ) also point out that 
they use a slightly different value of E{B — V) than in 
earlier studies, which could cause them to obtain a larger 
value of (m — M)q. 

We may also use our WD model calculations to mea- 
sure the extinction toward this cluster. The WD cool- 
ing models we have used provide us with expected col- 
ors for cluster WDs. These are compared with observed 
colors in Table [H Using only the three high-mass clus- 
ter member WDs, the weighted average E(B — V) = 
0.111 ± 0.025 and E(U - V) = 0.065 ± 0.028 for the 
cluster. If we include the lower mass potential cluster 
members, we obtain E(B - V) = 0.099 ± 0.021 and 
E(U - V) = 0.110 ± 0.025. Both valu es for E(B - V) 
are co nsistent with the value used by ISaraiedini et all 
(2004) (E(B - V) = 0.10 ± 0.01). The values we obtain 
for E(U — V) are lower than E(U — V) = 0.165 we would 
expect based on the E(B — V) value and our assumed 
reddening curve; however, given the uncertainty in our 
U -band calibrations (see Table [I} , we believe that this 
difference is not significant. 

5.4. The Initial-Final Mass Relation 

Having identified likely cluster member WDs, we now 
place them on the initial-final mass relation (IFMR, Fig- 
ure [T2"|) . The gray points show data from the litera- 
ture, with initial and final masses re-determined from 
published T e g and logg using our adopted evolution- 
ary models and methodology. References to the pub- 
lished data and the crucial assumptions are given in 
iWilliams et~"aT1 (|2008f) . The high- mass WDs selected for 
cluster membership (LAWDS 15, 17 and S2), represented 
by large black filled cir c les, f all directly on the trend 
found by iFerrario et al.1 (|2005l ) and represented by the 
dashed line. LAWDS 9 and 20 (intermediate-sized filled 
circles) fall well below the established IFMR, as does 
LAWDS SI (the candidate double degenerate). If cluster 
members, their relatively low masses compared to other 
WDs with similar progenitor masses would indicate that 
some mechanism, such as close binary evolution or en- 
hanced mass loss, is necessary to explain their location 
in the IFMR (see p~3]) . 

The three massive cluster mem ber WDs fall almost 
exactly on the linear IFMR from IFerrario et al.l (|2005l 
dashed line). For il lustration, we als o include the semi- 
empirical IFMR of Weidemann (200 (1 solid curve ) and 
the theoretical, Z = 0.019 IFMR of |Maiigd <|2001l . dot- 
ted curve). These relations bracket the empirical data. 
These points also begin to fill a region of the diagram 
(Mi n it = 4Mq — 5Mq) previously relatively unpopulated. 

The uncertainty in the age of the open cluster is a 
source of systematic error we now consider. As the life- 
time of the progenitor star is the cluster age minus the 
WD cooling age, an increase (decrease) in the assumed 
cluster age results in an increase (decrease) in the pro- 
genitor star lifetime, which corresponds to a decrease 
(increase) in the progenitor star mass for each WD. In 
other words, an older assumed cluster age leads to sys- 



tematically lower progenitor mass for every cluster WD, 
and vice versa. For this reason, we do not simply fold 
the cluster age uncertainty into the individual WD error 
bars, where it is likely to be misinterpreted as a random 
error for an individual star. 

The change in the derived IFMR due to the uncertainty 
in the age of NGC 1039 is shown in Figure [T2l as small 
points to the left and right of each cluster WD, which 
represent the IFMR for assumed ages of 250 Myr and 
200 Myr, respectively. For this cluster, the uncertainty 
in the star cluster age does not qualitatively affect the 
conclusion that the NGC 1039 IFMR agrees with the 
existing empirical IFMR. 

Figure Q2] shows the empirical IFMR for only those 
objects included in the LAWDS studies. There is con- 
siderably less scatter around the IFMR when only these 
points are plotted. This could be indicative of the het- 
erogeneity of the current open cluster WD sample, which 
is derived from several studies using a variety of instru- 
ments and atmospheric models. This suggests that fur- 
ther study into the inter-comparison of different samples 
is warranted. 

6. CONCLUSIONS 

We present photometric and spectroscopic data on the 
WDs in the field of the young open cluster NGC 1039. 
We use UBV imaging to photometrically select 44 objects 
as WD candidates. Of these 44, we spectroscopically 
identify 32 objects. Twelve are AGNs, one is an A-type 
star, 17 are DA or DAZ WDs, one is a DB WD (LAWDS 
26), and one is a DC WD (LAWDS 7). Five DAs are 
selected as possible cluster members- three high mass (> 
0.8-M©) WDs, and two low mass WDs. The membership 
of the DB WD remains in question and requires further 
analysis. The distance modulus of an additional DA is 
< 0.75 mag brighter than the distance modulus of the 
cluster, making it a potential (but not highly compelling) 
double-degenerate binary cluster member. 

The three high-mass cl uster members li e dire ctly on 
the established IFMR of IFerrario et "all (|2005h . while 
two lower-mass candidate members lie well below the 
relation. This could be due to binary evolution or 
metallicity-enhanced mass loss, though the evidence for 
either is weak. As the mass of these two latter WDs is 
similar to the mass at the peak of the field WD mass 
distribution, and as we expect to find ~ 2.4 field WDs 
meeting our membership criteria, these lower-mass WDs 
may not be true cluster members. Proper motion mea- 
surements will be necessary to confirm the membership 
of all five selected cluster member WDs. 

Based on the photometry of the five cluster mem- 
ber WDs, we obtain a cluster WD distance modulus of 
(m - M) v = 8.87 ± 0.10 and a cluster WD color ex- 
cess of E(B — V) = 0.10 ± 0.02. These values are consis- 
tent with recent measurements bv lSaraiedini et ail (|2004L 
(m - M) v = 8.98 with E(B-V) = 0.10). 

We note that the massive WDs at the distance modu- 
lus of the cluster all have cooling ages shorter than the 
age of the cluster. This indicates that results from main- 
sequence fitting do not seriously underestimate the clus- 
ter age. 

Eight out of 17 DA WDs have detectable Ca II K ab- 
sorption. The fraction of DAs showing Ca II absorption 
(47%) is significantly higher than the fraction that has 
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been found in other studies (~ 25%). Evidence suggests 
that the absorption is interstellar in origin, although this 
cannot be confirmed with present data. 

We use distance modulus and WD cooling time to se- 
lect probable WD cluster members, and while these cri- 
teria can provide us with a list of strong cluster member 
WD candidates, a study of proper motions of the cluster 
and these objects would significantly enhance our cer- 
tainty of their membership status. We leave this for a 
future work. 
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B-V 

Fig. 1. — B — V, V color-magnitude diagram for NGC 103 9. Left: our photometry alone; note that sources with V < 13 are saturated. 
Second left: our photometry plus that of cluster members from Jones & Prosscr ( 1996). Second right: combined photometry with isochrones 
and cooling curves overplotted using the [Jones & Prosscr (1996;) determination of the cluster (m — M )y. Right: combined photometry 
with isochrones and cooling curves adjusted to the cluster (m — M)y derived by Saraicdini ct al. (2004). In al l pane ls, small dots indicate 
our photometry and larger dots show high-probability proper motion cluster m embers from Jones & Prosscr (1996). The mam-sequence 
isochrones are a Z = 0.013, logT e ; = 8.35 P adova isochrone interpolated from lGirardi et al.l l|2002l , sol id), a Z = 0.013, log r c ; = 8.35 
Yonsei-Yale isochrone (Dcmarquc ct al. 2004, dashed), and the fiducial main sequence of NGC 2168 from lvon Hippel et al. (2002, dotted) 
adjusted to the metallicity of NGC 1039. Cooling curves for cluster WDs with cooling times less than 300 Myr and masses of O.4M0 (top), 
O.8M (middle) and 1.2M Q (bottom) are shown for DAs (solid) and DBs (dotted). 
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u-v 

Fig. 2. — U—V, V color-magnitude diagram for NGC 1039. Symb ols are as in Figure Q] Left: our photometry alone. Middle: our 
photometry with isochrones and cooling curves overplotted using the I Jones fc Prosscr (1996) determi nation of the clu s ter (m — M)y. 
Right: our photometry with isochrones and cooling curves adjusted to the cluster (m — M)v derived by Saraicdini ct al. ( 2004). 




Fig. 3. — B — V,U—B color-color diagram for NGC 1039. Cooling curves for cluster WDs with cooling times less than 300 Myr and 
masses of 1.2Mq (top), 0.8Mq (middle), and 0.4Mq (bottom) are shown for DAs (solid) and DBs (dotted). Vertical and horizontal dashed 
lines mark the boundaries of the WD candidate selection region. 




Fig. 4. — Lower left-hand regions of the color-magnitude diagrams for NGC 1039. Cooling curves are described in Figures [T] and [2] 
Dashed lines at the faint ends of the cooling curves mark the colors and V magnitudes of WDs which have been cooling for logr = 8.35 
(~ 225 Myr). Vertical and horizontal dashed lines mark boundaries of the WD candidate selection region. Open circles indicate objects 
that meet all of our WD selection criteria. 
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Fig. 5. — Upper left-hand region of the color-color diagram for NGC 1039. Cooling curves are described in Figure [3] Vertical and 
horizontal dashed lines mark boundaries of the WD candidate selection region. Open circles indicate objects that meet all of our WD 
selection criteria. 



14 



Rubin et al. 



8 



-LAWDS N3 



CD 
K 



2 - 



n — i — i — i — I — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 
LAWDS 7 



4*1^40101,^ LAWDS 18 



LAWDS 102 



I LAWDS SI ^ 



LAWDS 14 



LAWDS 19 






b=i i i i I i i i i I i i i i I i i — i — i i i i l i i i i l i i i i l i iz 

3500 4000 4500 5000 4000 4500 5000 

MA) MA) 



Fig. 6.— Spectra of WDs in the field of NGC 1039. The spectra of the DB (LAWDS 26) and DC (LAWDS 7) WDs are shown at the top 
of the left-hand panel, followed by the spectra of the DA WDs in order of decreasing WD T e g , from top to bottom and continuing through 
the right-hand panel. 
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Fig. 7. — Balmer-line profiles (light histograms) and best-fit models (heavy curves) for NGC 1039 DA WD spectra. The abscissa shows 
the offset in Angstroms from the center of each Balmer line. From bottom to top in each panel, the H/3 through H9 line regions are plotted. 
Each section of spectrum has been normalized to the pscudocontinuum just outside the plotted region and offset arbitrarily in the vertical 
direction. (A color version of this figure is available in the online journal.) 
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Fig. 8, — Histogram for all spectroscopically confirmed DA WDs in (m — M)y bins of width 0.3 mag. The I Jones & Prosser ( 1996) and 
Sarajcdini ct al. (200!) determinations of (m - M) v for NGC 1039 are marked with dashed lines. 
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(B-V) (U-V) 

Fig. 9. — WD candidate selection region of the color-magnitude diagrams for NGC 1039. Cooling curves are described in Figures [T] and 
[2] The symbols indicate high mass (> 0.8Mq) cluster member WDs (large filled circles), lower mass cluster member WDs (intermediate 
filled circles), the candidate double-degenerate WD (intermediate open circle), the DB (membership undetermined; open triangle), non- 
cluster member WDs (small filled triangles), non-WDs (crosses), unobserved WD candidates (small open circles), and stars not meeting 
our selection criteria (small points). Horizontal error bars show la uncertainties in de-reddened color, la uncertainties in My for WDs 
are smaller than the plotted point sizes. 
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(B-V) 



Fig. 10. — WD candidate selection region of the color-color diagram for NGC 1039. Cooling curves are described in Figure [3] symbols 
are as in Figure [9] Error bars show la uncertainties in de-reddened color. 
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Fig. 11. — Measurements presented inQylceicr ct al. (2005], grey triangles) with our Ca II K EWs and parallaxes (computed from our 
calculated (m — M)o values) overplotted. Our vertical error bars are smaller than the plotted point size. 2a EW upper limits on non- 
detections are also plotted. Our EWs fall within the scatter of the measurements of interstellar Ca II, sugg esting that Ca II in the WD 
spectra is interstellar i n origin. Dashed and dotted lines mark the parallax of NGC 1039 as determined by Sarajedini ct al. (2004) and 
IJones fc Prosscr ( 1996) respectively. 
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Fig. 12. — IFMR with new results over-plotted. Gray triangles show WDs included from the literature. High mass (> O.SMq) NGC 
1039 WDs are represented by large filled circles; low mass NGC 1039 WDs are shown with intermediate-sized filled circles. The open circle 
represents the candidate double-degenerate binary cluster WD. Error bars on NGC 1039 data points indicate the combined fitting errors 
for each point. Systematic changes in the location of the NGC 1039 WDs due to the uncertainty in the cluster ag e are s hown by the small 
points to the left (250 Myr) and right (20 Myr) of each WD. We include the linear IFMR fro m IFerrario et al.l l|2005l . dashed line), the 
semi-empirical IFMR of IWeidemanr] (|200d . solid curve), and the theoretical, Z = 0.019 IFMR of lMarigol l |2001l . dotted curve). 
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Fig. 13. — IFMR with only points from LAWDS plotted; symbols and lines are as in Figure [T2l The scatter appears to be smaller here, 
perhaps indicating that some of the scatter in the combined IFMR is due to the heterogeneity of the data set. 
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TABLE 1 

Transformation equation coefficients. 
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Adopted from Williams ct al. (2008) 



TABLE 2 

Photometry of candidate WDs in the field of KGC 1039. 
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TABLE 3 

Spectral fits for DA WDs in the field of NGC 1039. 
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Rubin et al. 



TABLE 4 

Progenitor masses for NGC 1039 WDs. 



LAWDS 


%ii a 


"Mvro a 


M init (M Q ) 


ID 


M 


M e 


log r cl = 8.30 log t c1 =8.35 log r ci = 8.40 



LAWDS 9 0.506 0.063 6.73 5 - 78 -a85 5 - 13 

LAWDS 15 0.870 0.074 4.72 4 - 41 ±o 22 4 - 14 

LAWDS 17 0.906 0.072 5.16 4 - 76 ±o.20 4 - 41 

LAWDS 20 0.561 0.072 18.38 8.39j;!j: % 6.51 

LAWDS SI b 0.578 0.070 4.37 4.14lg;^ 3.93 

LAWDS S2 0.844 0.074 4.23 4.021!! " 3.83 



a From Table [3p Potential binary member 



TABLE 5 

CA II K EW MEASUREMENTS FOR 
WDS IN THE FIELD OF NGC 1039. 



LAWDS 


Ca II K EW a dEW 


ID 


mA 


mA 


LAWDS 7 


< 147 




LAWDS 9 


280 


140 


LAWDS 14 


< 70 




LAWDS 15 


168 


70 


LAWDS 17 


< 300 




LAWDS 18 


240 


130 


LAWDS 19 


< 300 




LAWDS 20 


< 260 




LAWDS 22 


< 380 




LAWDS 25 


< 540 




LAWDS 26 


230 


100 


LAWDS 34 


320 


110 


LAWDS 40 


< 340 




LAWDS 102 


350 


100 


LAWDS N3 


180 


30 


LAWDS SI 


180 


50 


LAWDS S2 


160 


40 


LAWDS S3 


< 100 




LAWDS S5 


230 


110 


a Upper limits 


on EWs are 


at the 2(T 


level. 






TABLE 6 





Comparison of WD model colors with observed colors. 



LAWDS (m - M) v (B — V) c (B — V)o d E(B -V) (U - V) c (U - V) d E(U - V) 



LAWDS 9 a 


8.679 


0.119 


0.054 


0.065 


-0.416 


-0.648 


0.232 


LAWDS 15 


8.945 


-0.027 


-0.128 


0.101 


-1.094 


-1.210 


0.116 


LAWDS 17 


8.838 


0.078 


-0.108 


0.186 


-1.162 


-1.167 


0.005 


LAWDS 20 a 


8.931 


0.153 


0.084 


0.069 


-0.292 


-0.603 


0.311 


LAWDS SI b 


8.533 


-0.044 


-0.092 


0.048 


-0.919 


-1.061 


0.142 


LAWDS S2 


8.962 


-0.145 


-0.201 


0.056 


-1.328 


-1.374 


0.046 



a Potential low mass member b Potential binary member Observed color, from Table (2f Color determined 
from WD cooling models. 



